Available states and available space: Static 
properties that predict dynamics of confined fluids 
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Abstract. Although classical density functional theory provides reliable predictions 
for the static properties of simple equilibrium fluids under confinement, a theory of 
comparative accuracy for the transport coefficients has yet to emerge. Nonetheless, 
there is evidence that knowledge of how confinement modifies static behavior can 
aid in forecasting dynamics. Specifically, recent molecular simulation studies have 
shown that the relationship between excess entropy and self diffusivity of a bulk 
equilibrium fluid changes only modestly when the fluid is isothermally confined, 
indicating that knowledge of the former might allow semi-quantitative predictions of 
the latter. Do other static measures, such as those that characterize free or available 
volume, also strongly correlate with single-particle dynamics of confined fluids? Here, 
we investigate this question for both the single-component hard-sphere fluid and 
hard-sphere mixtures. Specifically, we use molecular simulations and fundamental 
measure theory to study these systems at approximately 10'^ equilibrium state points. 
We examine three different confining geometries (slit pore, square channel, and 
cylindrical pore) and the effects of particle packing fraction and particle-boundary 
interactions. Although average density fails to predict some key qualitative trends for 
the dynamics of confined fluids, we find that a new generalized measure of available 
volume for inhomogeneous fluids strongly correlates with the self diffusivity across 
a wide parameter space in these systems, approximately independent of the degree 
of confinement. An important consequence, which we demonstrate here, is that 
density functional theory predictions of this static property can be used together 
with knowledge of bulk fluid behavior to semi-quantitatively estimate the diffusion 
coefficient of confined fluids under equilibrium conditions. 
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1. Introduction 

Confining an equilibrium fiuid of particles to length scales on the order of several 
particle diameters changes both its static and dynamic properties. Classical density 
functional theory (DFT) can often make reliable predictions concerning the former, but 
implications of confinement for dynamics remain challenging to forecast for even the 
most basic models. For example, Enskog theory generalized for inhomogeneous fiuids [1] 
predicts that constant-activity confinement of the hard-sphere (HS) fiuid between 
hard walls will significantly decrease self-diffusivity parallel to the boundaries [2]. 
Recent molecular dynamics simulations of that system, however, show that this 
prediction is qualitatively incorrect, i.e., self diffusivity increases with constant-activity 
confinement [3j. Hydrodynamic theories, on the other hand, can predict how the 
presence of a single wall [1] or confinement between two walls [3 El [71 [8] impacts the 
self diffusion of a single Brownian particle in solvent. But it is not yet clear how 
to generalize these approaches to account for the effects that strongly inhomogeneous 
static structuring has on the transport coefficients of dense confined fiuids [Sj [9l [TO] . 
Given the absence of a reliable microscopic theory, the development of new qualitative 
or semi-quantitative heuristics for predicting dynamics of confined fiuids would be of 
considerable practical use. 

In this spirit, one productive strategy has been to first exhaustively simulate the 
static and dynamic behaviors of simple models of confined fiuids. The idea is that 
comprehensive study of these systems may reveal static quantities that strongly correlate 
with transport coefficients for a wide variety of confining environments. Knowledge 
of these correlations together with reliable predictions of the static quantities from 
equilibrium theory would, in turn, lead to indirect predictions for how confinement 
modifies dynamics. 

Indeed, recent simulation data covering hundreds of state points for various confined 
HS, Lennard- Jones, and square- well fiuids point to the existence of a robust, isothermal 
correlation between the self-diffusion coefficient D and the excess entropy per particle 
(over ideal gas) s*^^ — a relationship that is approximately independent of the degree of 
confinement for a wide range of equilibrium conditions [3l [Qj [Til [I2l [13]. The practical 
implication is that the mathematical form of the correlation for a given fiuid (obtained 
from bulk fiuid simulations) can be used together with knowledge of the excess entropy of 
the confined system (computed, e.g., via DFT) to make semi-quantitative predictions for 
how confinement will modify the self diffusivity. As has been discussed elsewhere [H [9] , 
this strategy can successfully predict subtle, confinement induced effects of packing 
frustration on relaxation processes, behavior not reflected in other static quantities 
traditionally thought to track dynamics, such as the average density. 

Although a fundamental and general derivation that explains the observed 
relationship between excess entropy and dynamics of conflned fluids is still lacking, the 
fact that the two are connected is not surprising. Excess entropy is a negative quantity 
that measures the extent to which static interparticle correlations reduce the number of 
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microstates available to the fluid. In fact, —s'^^ is often used as a metric for characterizing 
the "amount" of structural order present in condensed phase systems [T ^ [T5| [T6 | [T7 t [I8] . 
Since interparticle correlations strongly influence collisional processes, it makes intuitive 
sense that macrostate changes which increase structural order {—s'^^) might also tend 
to reduce single-particle mobility. Qualitative arguments such as these have previously 
been presented to rationalize empirically observed correlations between excess entropy 
and transport coefficients of both bulk [HI 1201 [21] and confined [21 El [TOl [IH [121 [IS] 
equilibrium fluids. Nonetheless, it is natural to wonder whether excess entropy is 
unique in this regard. Perhaps other measures also accurately forecast the implications 
of confinement for the dynamics of dense fluids. For example, does the mobility of 
inhomogeneous HS fluids also increase with the average amount of space available for 
particle motion? Do some measures of free or available volume correlate much more 
strongly with dynamics than others? Can one quantitatively, or semi-quantitatively, 
predict self diffusivity of confined fluids based on knowledge of how confinement affects 
the free or available volume of the system? 

Here, we put the above questions to stringent tests for a variety of confined 
fluid systems. To ascertain the effect of pore geometry on correlations between 
dynamics and thermodynamics, we study a monodisperse hard-sphere fluid confined 
to smooth hard-wall slit pores, square channels, and cylindrical pores. We also 
explore the effects of boundary interactions in the slit-pore geometry by examining 
cases for which square-well (attractive) or square-shoulder (repulsive) particle-boundary 
interactions are present. Finally, we investigate the dynamics and thermodynamics of 
a confined binary hard-sphere mixture which can be equilibrated in the fluid state at 
higher packing fractions without crystallizing than the corresponding single- component 
fluid. Together, this study represents, to our knowledge, the most comprehensive 
exploration of the relationships between static and dynamic properties in these basic 
inhomogeneous systems to date. We calculate the self-diffusion coefficient [via molecular 
dynamics simulations] , and excess entropy and various measures of available volume [via 
Transition Matrix Monte Carlo (TMMC) simulations and fundamental measure theory] 
at approximately 10^ state points. Our results illustrate that predictions of either excess 
entropy or a generalized measure of average available volume from classical density 
functional theory can be used together with knowledge of bulk fluid behavior to semi- 
quantitatively predict the diffusion coefficient of confined fluids across a wide range of 
equilibrium conditions. 

2. Connections between density, available volume, and dynamics: 
Preliminary evidence 

The available volume in a configuration of the single-component HS fluid refers to the 
space into which the center of an additional HS particle of equal size can be inserted 
without creating an overlap with existing particles or the boundary. It might also be 
considered a coarse measure of the space immediately available for particle motion in 
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that configuration of tlie system. For tlie bulk equilibrium HS fluid, the ensemble- 
averaged fraction of the total volume that is available (in the above sense) is given by 
Po = P/C [22], where p is number density, ^ = exp[/3/i]/A^ is activity, p is the chemical 
potential, f3 = [k-^T]'^, is the Boltzmann constant, T is the temperature, and A is the 
thermal deBroglie wavelength. It is known that increasing p monotonically decreases 
both pq and self diffusivity D for this system; i.e. dpo/dp < and dD /dp < (see, e.g. 
[18]), which is consistent with the intuitive idea that these static and dynamic properties 
might be linked. 
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Figure 1. Schematic of an iniiomogeneous HS fluid confined between boundaries in a 
slit pore geometry. The walls are placed a distance H = h + a apart, where h is the 
length accessible to particle centers. Dark regions indicate hard spheres and confining 
walls. Additional particle centers are excluded from the grey (overlap) region. The 
white region indicates the volume available for inserting an identical hard sphere of 
diameter a. 



Are density, available space, and self diffusion connected in the same qualitative way 
for inhomogeneous HS fluids? Previous studies have provided some information useful 
for addressing this question. For example, it is known that confining the equilibrium HS 
fluid between hard walls (while maintaining fixed ^) significantly increases the average 
particle density, i.e. (dph/dh^^)^ > 0, over a wide range of ph and /i [21 |3]. Here, 
Ph = N/ {Ah) is average density, is the total number of particles, A is the area 
of the interface between the fluid and one hard wall, and h is the center-accessible 
width of the slit pore [i.e., not including the "dead space" that the particle centers 
are excluded from due to their interaction with the boundaries (see Figured])]. The 
fact that Ph increases upon constant-,^ confinement initially appeared consistent with 
earlier kinetic theory predictions that D of this system would show a corresponding 
decrease [2] . However, recent simulation data have demonstrated that the kinetic theory 
predictions were qualitatively incorrect [3]. That is, both ph and D typically increase 
upon constant-,^ confinement [{dD / dph)^ > 0, the fluid gets "denser" on average and 
particles diffuse more rapidly], a trend that is the opposite of what might be expected 
based on the bulk HS fluid behavior. Interestingly, confined HS fluids show a different 
trend when an alternative thermodynamic constraint is applied. Specifically, if h is held 
fixed, then increasing ph has the effect of decreasing D, i.e., {dD / dph)h < 0] [3],[II]. This 
preliminary data illustrates that knowledge of how ph changes is not, in and of itself, 
enough to even qualitatively predict the implications of confinement for the dynamics of 
a fluid. Indeed, in Section Hlof this article, we present extensive numerical evidence for a 
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variety of confined fluid systems which underscores this point. We also explore whether 
adopting a definition of average density based on the total rather than center-accessible 
volume of the pore (see also [31 [H]) improves predictions for how confinement modifies 
dynamics. 

Does available volume show a more reliable correlation to dynamics than average 
density? The fractional available volume in an inhomogeneous fiuid is inherently a 
local quantity, and it is given by Poi{z) = pj(z) exp[/3uj"(z)]/,^j, where Pi{z) and u^{z) 
represent the singlet (one-particle) density and the wall-particle interaction potential for 
species i, respectively, evaluated at a distance z from one wall [23l[2ll[25]. The volume 
averaged quantity can be expressed as 



where the integral is over the particle-center- accessible volume Vc^i. For the special case 
of a single-component HS fiuid confined between smooth hard walls, Eq. ([T]) reduces to 
Po = Ph/C Note that since {dpo/dh~^)^ = {dph/dh~^)^, and {dph/dh~^)^ > across 
a wide range of ph and /i [21 [3], it follows that {dpo/dh~^)^ > for those conditions. 
This increase in the fraction of available space with constant-.^ confinement provides a 
simple physical explanation for the counterintuitive observation that D increases along 
the same thermodynamic path. The density-dependent behavior of the po under the 
constraint of constant h is also qualitatively consistent with the dynamical trends of the 
confined HS fiuid. In particular, both {dpo/dph)j^ < and {dD / dph)^ < 0. All of this 
strongly suggests that po is a more relevant static metric for particle mobility than the 
average particle density ph- 

In the following sections, we test the generality of these preliminary observations 
by carrying out an extensive quantitative comparison of the correlations between 
self diffusivity D and various static measures (density, excess entropy, and fractional 
available space) for single-component and binary HS fiuids confined to a variety of 
channels with different geometries and particle-boundary interactions. The results 
clarify which of these static quantities reliably predict the implications of confinement 
for single-particle dynamics. 

3. Methods 

We study single-component HS fiuids of particles with diameter a both in the bulk 
and confined to channels with three types of geometries: (i) quasi-two-dimensional slit 
pores, (ii) quasi-one- dimensional square channels, and (iii) cylindrical pores. Specifically, 
we consider (i) seven slit pores with thickness H/a = 5,6,7,8,9,10, and 15 in the 
confining z direction (see Figure ^ together with periodic boundary conditions in 
the X and y directions, (ii) seven square channels with total x — y cross-sectional 
dimensions of H"^/ = 25, 36, 49, 64, 81, 100, and 225 together with a periodic boundary 
condition in the z direction, and (iii) six cylindrical channels of total diameter H/a = 
6, 7, 8, 9, 10, and 15 together with a periodic boundary condition in the axial z direction. 
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We take the interaction of a particle with a channel wall u^{s) in all cases to have 
a generic square-well form: 

( oo s < a 12 



where s is the shortest distance between the particle center and the wall of interest. 
For all three geometries, we study the case e„ = 0, i.e. smooth hard boundaries. 
Additionally, for the slit pore with size H/a = 5, we investigate cases with ew = 2kBT 
(a repulsive shoulder) and e„ = —2k-QT (an attractive well). 

We also consider a binary HS mixture confined between smooth hard walls in a 
slit pore of width H/ai = 5. For this system, the particle diameter ratio is taken 
to be o"2/o"i = 1.3 and the particle masses are proportional to their volume, i.e., 
1712 /mi = (o"2/cri)^- These parameter values closely mimic those examined in recent 
experiments of binary colloidal mixtures under confinement [26] . 

To explore dynamic properties of these fluids, we perform molecular dynamics 
simulations using an event-driven algorithm [27] in the microcanonical ensemble with 
N = 4000 particles for monodisperse HS systems and = 3200 for the binary HS 
systems. For bulk systems, we use a cubic simulation cell of volume V. For the confined 
systems, we adopt a rectangular parallelepiped simulation cell of volume V = H^HyHz 
with appropriate boundary conditions depending on geometry. We extract the self 
diffusivity D by fitting the long-time {t » al^Jm^) mean-squared displacement to 
the Einstein relation (Ar^^) = 2dpDt, where (Ar^^) corresponds to motions in the 
dp periodic directions. For the sake of clarity, we reserve the symbol D for the self 
diffusivity of fluids under confinement and -Dbuik for the self diffusivity of the bulk 
fluid. To obtain reliable estimates, we average self diffusivities over four independent 
trajectories. For simplification, we report quantities from this point forward implicitly 
non-dimensionalized by appropriate combinations of the characteristic length scale, 
taken to be the HS diameter of the smallest particle in the fluid cxi, and a characteristic 
time scale, given by aiy/miP. Thus, all energies are implicitly per unit k-oT, and T is 
effectively scaled out of the problem altogether. 

We obtain thermodynamic properties using grand-canonical transition-matrix 
Monte Carlo (GC-TMMC) simulation. For pure fluids, we use an algorithm presented 
by Errington [28] and for binary mixtures we employ a strategy developed by Shen 
and Errington p9], wherein one combines a series of semigrand ensemble simulations 
to construct the system's free energy over a wide range of densities and compositions. 
We conduct GC-TMMC simulations in a standard grand canonical ensemble where the 
volume V, temperature T, and activities {^5^2} are held constant and the particle 
numbers {Ni,N2} and energy E fluctuate. For notational convenience, we denote the 
sets {A^i,iV2} and {^5^2} as N and ^, respectively, using conventional vector notation. 
The activity of component i is defined as = A^"^ exp(/ij), where /ij is the chemical 
potential and Aj is the thermal de Broglie wavelength. For the pure component GC- 
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TMMC simulations we present here, we set ^ = 1 and adjust the particle center- 
accessible volume Vc to make the total volume V ~ 1000. Simulations of the binary 
mixture use V = 125 and V = 245 for the bulk and confined fluids, respectively. For 
the bulk and confined simulations, we set = 173.7 and ^2 = 381.5. 

The key quantity we extract from the GC-TMMC simulations is the particle 
number probability distribution n(N). Once we obtain this distribution, we use 
basic statistical mechanics principles and histogram reweighting [301 to evaluate 
thermophysical properties over a range of activity values. First, we use histogram 
reweighting to deduce n(N) at a set of activities ^^ew generally different from that 
of the GC-TMMC simulation ^^j^, 

lnn(N; = lnn(N; ^si^) + J] iV,(ln6,new - In^^.sim), (3) 

i 

where it is understood that the probability distributions are not normalized. We obtain 
average particle numbers (N) from first-order moments of n(N), 

(N(0) = 5^Nn(N;0/5^n(N;0. (4) 

We calculate p and ph via normalization of (N) by V and K, respectively. 

We define excess entropy as the difference between the fluid's entropy and that of 
an ideal gas with the same density profile. The particle number probability distribution 
provides the density and composition dependence of the Helmholtz free energy at a given 
temperature. Therefore, we combine knowledge of n(N), average excess configurational 
energies If^^CN), and particle-number-specific spatial density distributions p(N, r) to 
obtain the total excess entropy S"^^ [TH [131 IB] , 

5'=^(N) = f/'=^(N) + lnn(N)/n(0) 

+ ^|lniVJ-iVanei-iVaniVi + ^ pi(N, r) lnpi(N, r)dr| (5) 

We also predict the thermodynamic quantities of the bulk single-component and 
binary HS fluids using the Carnahan-Starling [31j and Boublik-Mansoori-Carnahan- 
Starling-Leland [32l [33] equations of state, respectively. We predict the thermodynamic 
properties of confined HS fluids using a recent modification [M] of Rosenfeld's 
fundamental measure theory [35]. Fundamental measure theory is a classical DFT 
of inhomogeneous fluids that has been shown to accurately predict structure and 
thermodynamics of confined HS systems in various restrictive geometries up to very 
high densities [36]. For numerical evaluation of the DFT for slit and cylindrical pores, 
we use Picard iterations on a grid spacing of 0.005. We update densities according to, 
[p]n+i = 0-95[p]^ + 0.05[p]°"*, where [p]n+i and [p]™ are the input density profiles at 
the n + 1^^ and n*^ iterations, respectively, and [p]^"* is the output density profile at 
the n*^ iteration. We stop Picard iterations when the relative change in output density 
profile between two successive cycles [1 cycle = 20 iterations] becomes less than 10~^. 
For numerical evaluation of the DFT in the square channel geometry, we use Sandia 
National Laboratories Tramonto package [37] . We adopt a grid of 0.05 x 0.05 for H = 5,6 
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and a grid of 0.075 x 0.075 for H > Q. We stop the minimization algorithm when the 
relative or absolute change in the grand potential is less than 10~^. 

4. Testing structure-property relations for predicting self diffusivity of 
confined fluids 

In this section, we explore the accuracy of the following strategy for predicting the self 
diffusivity of confined HS fluids: (1) determine the value of a static quantity x of a 
confined fluid believed to be relevant for dynamics (e.g., its density, excess entropy, or 
fractional available volume), and (2) input this value into the relationship between self 
diffusivity and that same static quantity for the bulk fluid, D^^^^^x), to estimate the 
confined fluid self diffusivity, D. Of course, such a strategy can only provide approximate 
predictions. While there is a one-to-one relationship for the equilibrium HS fluid between 
the self diffusivity and any one of the aforementioned static quantitiesjj] the dynamic 
properties of the confined fluid generally depend on a larger number of variables (e.g., the 
dimensions of the confining geometry, the nature of the particle-boundary interactions, 
and the chemical potential). Nonetheless, the hope is that one can discover a static 
quantity x whose relationship with D is largely insensitive to the effect of confinement. If 
so, the bulk structure-property relationship D^^^^^ can be used to predict D independent 
of the other parameters of the confined system. Systematic tests of this idea should 
give new insights into the structural properties that are most relevant for single-particle 
dynamics of inhomogeneous fluids. 

To investigate the accuracy of predictions by this approach, we use the "exact" 
results of molecular simulations to examine the ratio of the bulk fluid self diffusivity 
to that of confined fluids with the same value of x, i.e., D^^^^k/-^- Since one is often 
interested in both D and the effective characteristic time associated with diffusive motion 
{D~^), we present plots of D^^y^/ D in this work on a logarithmic scale, a representation 
for which overpredictions and underpredictions of D by the same factor appear the same 
distance from unity. We also present statistics associated with the relative errors of the 
predictions for different x. For each x that we consider here, we analyze roughly 10^ 
state points of the equilibrium HS fluid confined to the various pore geometries described 
in Section [31 This data, when taken as whole, spans approximately four decades in D. 

4.1. Confinement in channels with smooth hard boundaries 

Here, we examine the ratio D^^y^/ D for the single- component HS fluid confined to 
various geometries by smooth hard boundaries [i.e. = 0, see ([2])]. We begin by 
testing the predictions that follow from assuming that x = ph, the number density 
averaged over the particle-center-accessible volume of the pore, is the relevant static 
metric for dynamics [see Figure [2]^ a-c)]. It is immediately clear from the data that ph 

I This is true if the self diffusivity is appropriately non-dimcnsionalized, as we have done here, to 
remove the trivial effect of the thermal velocity of the particles. 
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Figure 2. Ratio of self diffusivity of a bulk HS fluid to that of a confined HS fluid with 
the same value of a static quantity x, D^^^k/D, obtained via molecular simulations. 
Data are shown for the fluid confined to slit pore [a, d, g, j], square channel [b, e, h, 
k], and cylindrical pore [c, f, i, 1] geometries. The static quantity x is indicated in 
each of the four panels. 20% bounds on relative error in confined fluid self diffusivity 
"predictions" by using the bulk structure-property relation D^^^y. are shown by red 
dash-dotted lines. The equilibrium fluid states shown here span the density range 
< Ph < po, where po ~ 0.95 for (a-c) and 1 < po < 1-25 for (d-1), depending on 
pore size. Pores shown have confining dimensions of 5 (O ), 6 (□), 7 (0), 8 9 (x), 
10 (*), and 15 (a). These dimensions correspond to channel width H for slit pores 
and square channels and channel diameter d for cylindrical pores. All channels have 
smooth hard boundaries. 



does not, in itself, provide a good basis for prediction. HS fluids confined in slit-pore, 
square-channel, and cylindrical geometries generally exhibit a wide range of D for each 
Ph, with the fastest dynamics occurring in the smallest pores. In fact, note that the bulk 
structure-property relation D^'^^y_ can underpredict D by nearly a factor of ten for fluids 
in the most restrictive geometries. The performance of the bulk structure-property 
relation using x = ph is actually even worse than it appears in Figure [2](a-c) for the 
following reason. The freezing transition occurs at a density of 0.945 for the bulk HS 
fluid, which provides an upper limit on values of ph that can be used for predictions 
using ^buijj. However, center-accessible densities for the equilibrium fluid in the smallest 
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square- channel and cylindrical pores can reach as high as ph ~ 1.25. Thus, the bulk 
structure-property relation D^'^^^ cannot even make predictions for a significant fraction 
of the equilibrium state points for highly confined HS fiuids. All of this confirms the 
preliminary expectation discussed in Section [2] that knowledge of ph and bulk fiuid 
behavior is not enough to predict the self diffusivity of confined fiuids. This conclusion 
is consistent with the earlier observations of Mittal et al. [11] concerning a smaller set 
of data for the HS fiuid confined to slit pores. 

In Figure E] (d,g, and, j), we again present D^^^^/ D for the HS fiuid confined to 
slit pores, but now -D^uik corresponding bulk fiuid relation between diffusivity 

and one of three alternative static properties (x): average density based on total pore 
volume p = pfi{l — H~^), excess entropy per particle s"^^, and fraction of available volume 
Pq. The data in these plots corresponds to confined fiuid states with packing fractions 
that vary from the dilute gas to the freezing transition for pore widths H > 5. As 
can be seen, each of these static measures can provide semi-quantitative predictions for 
confined fiuid diffusivities when used together with the corresponding bulk structure- 
property correlation. In fact, for 93% {x = p), 97% {x = s"^^), and 100% (x = po) of 
equilibrium state points for these systems, the predictions provided by -D^uik are within 
20% of the "exact" MD data for D. Note that the very small fraction of overpredictions 
based on p or s'^^ that exceed 20% relative error correspond to the high density, low D 
state points near the freezing transition. 

Based on the slit-pore data, it might be tempting to conclude that total-volume- 
based average density p tracks dynamics nearly as reliably as s*^^ and po for confined 
fiuids. To provide a more stringent test of this preliminary conclusion, we now examine 
-^buik/-^ for HS fiuids confined to quasi-one-dimensional square channel and cylindrical 
pore geometries with edge dimensions H > 5 and diameters d > 6, respectively. Fluids 
confined in these geometries have a significantly higher percentage of particles near the 
boundaries than in the corresponding slit pores, and hence the effects of confinement 
on both structure and dynamics should be more pronounced. 

Figure [2](e,f) shows that for square channel and cylindrical geometries, self- 
diffusivity predictions based on D^^^y_ can be significantly higher than the actual D 
of a confined fiuid with the same p. In fact, for the densest fiuid systems studied here, 
the bulk structure-property relation D^^^y_ is between 2 and 4 times larger than D, 
depending on H. Furthermore, Figure [2]^h,i,k,l) illustrates that D^^^^^ predictions for 
X = s'^^ or X = po are in the semi-quantitative range for a larger fraction of state points 
than those based on x = p. Specifically, -Dbuik within 20% of D for 46% (x = p), 82% 
{x = s^^), and 95% {x = po) of the state points. The main differences occur for high 
density, low D state points, where predictions based on fractional available volume are 
significantly more accurate than those based on excess entropy or density. 

Another relevant test case for comparing which of p, s'^^, or po most faithfully tracks 
dynamics is to vary the degree of confinement while fixing p, an idea motivated by an 
earlier study by Mittal et al [3] . In particular, Mittal et al demonstrated that D and s^^ 
of a HS fiuid oscillate in phase when H of the confining slit pore is varied (for H < 5) and 
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Figure 3. Self difFusivity D of a HS fluid confined in narrow slit pores of width 
iJ = 2 — 5 by smooth hard boundaries. The density of the confined fluid is fixed at 
p — (6/7r)0.4. We compare molecular dynamics simulation data (O ) for the confined 
fluid with the self diffusivity of a bulk HS fluid, -D^uik' same value oi x = p 

(green dash-dotted line), x = (red square), and x = pq (blue triangle). Solid lines 
are shown as a guide to the eye. D and s"'^, calculated via molecular dynamics and 
TMMC simulations, respectively, are taken from Fig. 7 of 3 . po was calculated via 
TMMC simulations described in Section [3] of the present work. 



p is held constant. The maxima in D (high particle mobility) and s'^^ (weak interparticle 
correlations) occur for integer values of H, geometries which naturally accommodate the 
layering of particles near the boundaries. The minima in D (low particle mobility) and 
s'^^ (strong interparticle correlations) occur for non-integer values of H, which frustrate 
this natural layering pattern. Along similar lines, Goel et al [9] recently demonstrated 
that particle-boundary interactions that flatten the density profile of a confined fluid 
generally reduce D and s*^^, while those which increase layering can have the opposite 
effect. All of this suggests that excess entropy captures some of the subtle frustration 
induced effects that layering has on both interparticle correlations and single-particle 
dynamics [3]. Does po also capture these effects? A very recent study of Mittal et al [TO] 
suggests that it might. In particular, the authors of that study showed that the local 
fraction of available volume Pq{z) and the position-dependent self diffusivity normal to 
the boundaries Z)_l {z) of a confined HS fluid are highest in regions of high local density 



Figure [3] provides a more direct test of this idea. In particular, it shows the D data 
of Mittal et al [3J for a HS fluid confined between hard walls calculated via molecular 
dynamics simulations. We have also included on this plot predictions from the three 
bulk structure-property relations -Dbuik; where x = p, x = s''"', and x = po- Since p is 
fixed here, it is evident that D'^^^^ is not able to predict the oscillatory trends in the 
dynamics data. However, note that both D^^^^^ and -Dbuik pi'edict the correct oscillatory 
behavior. In fact, the predictions of -0^°!^ virtually quantitative over the entire range 



p{z). 
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0.001 



Figure 4. Ratio of self diffusivity of (a) small particles and (b) large particles of a bulk 
binary HS fluid mixture to that of a corresponding confined HS fluid mixture with the 
same value of a static quantity x, D\^^\]^/ obtained via molecular simulation. The 
static quantity used for making predictions isa; = p(0),a:: = (□), and x = po^ 
(a). Data are shown for the H = 5 slit-pore geometry with the total packing fraction 
in the range 0.025 — 0.52. The mole fraction of the small spheres is 0.75. 20% bounds 
on relative error in self-diffusivity predictions are shown by red dash-dotted line. 

Should we expect /^^"j^ to generally track the dynamics of dense, confined HS 
fluids more accurately than i^bu'ik'^ other words, what is more relevant for dynamics 
of inhomogeneous fluids: available space or available states? We further explore that 
question by examining the behavior of the confined HS mixture discussed in Section |3l 
Specifically, by studying this binary fluid mixture in a slit pore with if = 5, we are 
able to probe confined fluid states with packing fraction as high as 0.52 (compared 
with the highest packing fraction of 0.46 for a monodisperse HS fluid confined in a 
slit-pore of ii = 5) At = 0.52, the fluid already exhibits dynamic signatures of 
supercooling, e.g., the emergence of a plateau in the time dependence of the mean- 
squared displacement. The corresponding D at this packing fraction (= 0.002) is an 
order of magnitude smaller than the smallest D for the confined monodisperse fluid 
(= 0.02) in the slit pore geometry. 

Figure IH^a,b) shows the ratio of the bulk self diffusivity to confined self diffusivity, 
-^buik/-^) for the small and large particles of the mixture, respectively. Again, the 
comparisons are made to the bulk fluid mixture of the same composition and density 
[x = p), excess entropy {x = s^^), or fractional available volume of the corresponding 
species {x = Poi). Note that the bulk structure-property predictions for self diffusivities 
of small and large particles are semi-quantitative (within 20%) for i^ > 0.1 based on 
any of the three aforementioned static quantities. However, relative errors in predictions 
based on p or s*^^ begin to increase sharply for D < 0.1. On the other hand, predictions 
based on p^^ remain semi-quantitative for all D > 0.02 (covering three decades in 
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D), with significant overpredictions occuring only for tlie densest three state points 
investigated. Thus, it appears that, for single-particle dynamics, fractional available 
volume is the most relevant of the three static measures investigated here. The question 
of whether there exists an alternative static measure x, such that -D^^jj^ tracks D for 
deeply supercooled mixtures, is currently an open one. The answer to that question will 
likely have important implications for understanding how confinement shifts the glass 
transition of fluids. 



4-2. Particle-boundary interactions and a generalized measure of available volume 



Thus far, we have only considered the geometric (i.e., packing) consequences of 
confinement on dynamics. How does the presence of finite particle-boundary interactions 
affect the picture described in the previous section? We explore the answer to this 
question by studying a monodisperse HS fluid conflned to a slit pore geometry by 
smooth walls with either square-shoulder (repulsive) or square-well (attractive) particle- 
boundary interactions [for details, see Eq. 
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Figure 5. Ratio of self diffusivity of a bulk HS fluid to that of a confined HS 
fluid with the same value of a static quantity x, D^^^^/D, obtained via molecular 
simulation. Data include systems with (a) square-shoulder (repulsive) and (b) square- 
well (attractive) particle-boundary interactions (see Eq. [J) with H = 5 in the slit pore 
geometry. The static quantity used for making predictions is x — p {O )^ x = s'^^ 
(□), X ^ po (a), and x — pow (A). 20% bounds on relative error in self-diffusivity 
predictions are shown by the red dash-dotted line. Particle center accessible density for 
the fluid spans the range < ph < 1.05. po overpredicts high-diffusivity state points 
in panel b by more than 1000% and those data points are off the scale of the graph. 



Figure [5](a and b) shows the ratio of the bulk to confined self diffusivity, D^^^^^/D, 
for these two cases, respectively, with comparisons between bulk and confined fiuids 
being made at the same density (x = p), excess entropy (x = s'^^), and fractional 
available volume (x = po) . Clearly, the fractional available volume fails to track the 
dynamics of the confined fiuid for both cases presented. However, this should perhaps 
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be expected. Since the interactions of the particles with the boundaries in these cases 
are strongly position-dependent, all free space is not equally "available" to the particles. 
To account for this energetic imbalance, we suggest a generalized available volume (pow) 
that appropriately weighs the local available space by the Boltzmann factor of the 
particle-boundary interaction. 



where, in all cases, the superscript "ig" denotes the corresponding quantity for an ideal 
gas confined to an identical slit pore. As one can see from the above equation, since this 
generalized available volume inherently relates the thermodynamic state of the confined 
fluid to that of an ideal gas, it bears some resemblence to an excess thermodynamic 
property. 

It is important to point out that this generalized available volume has several 
distinguishing features, (i) It reduces to Pq in the limit of fluids confined to smooth 
hard boundaries. Thus, all of the results presented earlier in this paper for po will 
remain unchanged for those systems if one instead uses pow- (h) Unlike po or density, 
there is no arbitrary choice that needs to be made about the volume over which one 
should do the averaging. This quantity is the same no matter whether averaging is 
carried out over the center-accessible or the total volume of the fluid. Put differently, 
this definition removes any arbitrariness about the effective "diameter" of the fluid- 
boundary interaction, (iii) The quantity po^ can be computed directly from knowledge 
of standard thermodynamic and system parameters, namely p, ^, and Uv,{z). (iv) It is 
not limited to HS fluid systems. In fact, the computation of po^ from Eq. (E]) does not 
even require information about the particle-particle interactions, as long as the other 
thermodynamic quantities can be measured. 

How well does this new generalized measure of available volume track self diffusivity 
when finite particle-boundary interactions are present? Figure [5] clearly illustrates that 
p^ corrects for the problems that p^ faces in predicting the dynamics in these cases. 
In particular, the maximum error in self-diffusivity predictions based on p^ is 16% 
across the entire range of packing fractions investigated, which makes it a more reliable 
predictor of single-particle dynamics than either s*^^ or p for these systems. As was seen 
earlier, diffusivity predictions based on the bulk structure-property relation for pow are 
considerably more accurate than that for p or s*'^ when considering high density, low 
diffusivity state points {D < 0.1). 

5. Using DFT together with bulk structure-property relations to predict 
dynamics of confined fluids 

Above we have shown that knowledge of -Dbui^ for the bulk HS fluid (x = s"^^ or x = po^) 
together with the value of x in confinement is enough for semi-quantitative prediction 
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Figure 6. Ratio of self diffusivity of a bulk HS fluid to that of a confined HS fluid 
with the same value of a static quantity x, D^^^^^/ D. Here we show a comparison of 
^buik/^ with X obtained from the DFT (lines) or TMMC (symbols) calculations, (a, 
b, and c) Data for a monodispcrse HS fluid confined by smooth hard boundaries to a 
slit pore [H — b], a cylindrical pore [iJ = 6] , and a square channel [H — 5] geometry, 
respectively, (d and e) Data for small and large particles of a binary HS mixture 
confined in a iJ = 5 slit pore, respectively, (f and g) Data for a monodispcrse HS 
fluid confined between the repulsive and attractive boundaries discussed in the text, 
respectively, placed at a separation of = 5 in slit pore geometry. The density range 
for the state points shown are the same as in Figure [5] for [a-c], Figure S] for [d and e], 
and Figure [5] for [f and g] . 20% bounds on relative error in self-diffusivity predictions 
are shown by the red dash-dotted line. 



of confined (equilibrium) fluid diffusivity, D, across a wide range of parameter space 
in these systems. Thus far, we have used TMMC simulations to determine x for each 
confined fluid of interest. That raises the following question. Is classical DFT accurate 
enough in its predictions of x that one can eliminate the step of simulating the confined 
fluid altogether? In this section, we take a first step toward addressing this question. 
In particular, we present calculations of the ratio of self diffusivity of a bulk HS fluid 
to that of a confined HS fluid with the same value of x, D^^^^/ D, where x = s"^^ or p^. 
In each of the confined fluid is obtained directly from predictions of Rosenfeld's 

fundamental measure theory [3lj, an accurate DFT for these systems. 

Figure [6] shows relative errors in self-diffusivity predictions based on s^^ and pow, 
comparing cases with knowledge of the "exact" value of x in confinement (calculated via 
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TMMC simulations) and the predicted value of x (calculated via DFT). Selected cases 
explored in the previous section involving the three confining geometries (slit pore, 
cylindrical pore, and square channel), small and large particles of binary mixtures, 
and finite fiuid-boundary interactions are presented. In all cases, the D^uik/-^ curves 
obtained via the two routes (TMMC versus DFT) are virtually indistinguishable over 
the entire density range of equilbrium fiuid, a demonstration of the reliability of DFT 
for computing the static properties of these systems. As a result, it is clear that one 
can use the bulk structure-property relations discussed above together with predictions 
of X from DFT to make semi-quantitative estimates of confined fiuid self diffusivity for 
a wide variety of hard sphere systems. 

6. Conclusions 

Fluids trapped in small spaces feature prominently in science and technology, and 
understanding their properties is key for a number of research areas that range from the 
design of membranes to the engineering of microfiuidic devices. The static and dynamic 
properties of these confined fiuids can be very different than those of bulk samples. While 
quantitatively accurate theories like DFT are available for predicting static properties 
of confined fiuids, making even qualitative predictions for dynamics of inhomogeneous 
fiuids has long been a challenging endeavor. In this paper, we demonstrate how 
semi-quantitative (albeit indirect) predictions of dynamics from first principles are still 
possible, even in the absence of a theory, once one recognizes that certain relationships 
between static and dynamic properties are insensitive to confinement. 

This study provides a systematic and quantitative investigation of such 
relationships. In particular, we present a comprehensive study of the effects of 
confinement on the correlation between self diffusivity (D) and various thermodynamic 
measures for confined HS fiuids: particle center-accessible-volume-based and total- 
volume-based average densities {ph and p, respectively), excess entropy (s^""), and two 
average measures of fractional available volume (po Po^)- Our main findings are 
as follows. The bulk structure-property correlation, D^^^^, based on the first density 
measure, x = ph, severely underestimates D when ph of the confined fiuid is used as 
the input. Further, for dense confined fiuids, ph is often larger than the freezing density 
of the bulk fiuid, eliminating altogether the possibility of using the corresponding bulk 
structure-property relation for predictions. Self-diffusivity predictions based on the 
relation with total-volume-based density x = p provide a significant improvement over 
those involving ph, substantiating the earlier idea that p might be considered as a more 
natural measure of density for predicting dynamics [TT] . 

However, when one considers a wider variety of geometries [slit pore, cylindrical, 
square channel], confined fiuid mixtures, finite particle-boundary interactions, and a 
wide range of packing fractions, one finds that structure-property relations -D^ui^ based 
on excess entropy {x = s^^) and a new generalized measure of available volume {x = pow) 
provide much more accurate estimates for D than those based on either of ph or p. 
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Further, self-diffusivity predictions based on are significantly more accurate than 
those based on s"^^ under conditions of high packing fractions (e.g., supercooled fluids) 
and highly restrictive confining geometries (e.g., quasi-one-dimensional channels). The 
generalized available volume po^ may also be easier to compute based on experimental 
quantities than s*^^, since the former is related to average density and chemical potential 
in a simple way (see Eq. [6]). Importantly, neither p^ nor s^^ require one to arbitrarily 
define an averaging volume for the confined system. 

Predictions of x via classical DFT are accurate enough for inhomogeneous HS 
fluids that one can use them, together with the bulk structure-property relation 
(-^buik)' to make semi-quantitative estimates of confined fluid diffusivities. This 
effectively eliminates the need for simulating the confined fluid altogether, which 
might be particularly convenient in applications where one needs to estimate the 
dynamics of systems across a wide array of parameter space. For example, in the 
design of microfluidic systems, one might hope to screen a large range of possible 
particle-boundary interactions or confining geometries against design considerations. 
A preliminary application of this idea [9j is to use DFT to passively tune the transport 
properties of a confined fluid, in a controlled way, by modifying the geometry or 
boundary-particle interactions of the confined space. 

Can the aforementioned static measures predict the dynamics of fluids with 
continuous intermolecular potentials and/or attractive interactions? As mentioned 
earlier, recent data from molecular simulations have shown that there is indeed an 
isothermal correlation between the self-diffusion coefficient D and the excess entropy 
for a variety of confined fluids (e.g., Lennard- Jones, square- well. Weeks- Chandler- 
Andersen), approximately independent of the degree of confinement for a wide range 
of equilibrium conditions [121 E] ■ We are currently exploring the viability of using the 
generalized available volume for predicting dynamics in these fluids, and we will report 
on our findings in a future study. 
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